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We show that a recently proposed derivation of Bose-Einstein correlations (BEC) by means 
of thermal Quantum Field Theory, supplemented by operator-field evolution of Langevin type, 
allows a deeper understanding of a possible coherent behaviour of the emitting source and a clear 
identification of the origin of the observed shape of the 2-particle BEC function C2(Q)- We explain 
the origin of the dependence of the measured correlation radius on the hadron mass. The lower 
bound of the particle source size is estimated. 
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In our previous papers [1,2] we have established the master equation of the field-operator evolution (Langevin- like 
CN ' [3]) equation which allows one to gain a better understanding of possible coherent behaviour of the emitting source 
, of elementary particles. A clear identification of the shapes of both two-particle both Bose-Einstein and Fermi-Dirac 
correlation functions was observed in LEP experiments ALEPH [4], DELPHI [5] and OPAL [6] which also suggested 
pLn [ a dependence of the measured correlation radius on the hadron mass m. 

We focused [2,1] our attention on two specific features of Bose-Einstein correlations (BEC) clearly visible when BEC 
| are presented in the language of Quantum Field Theory (QFT) supplemented by the operator-field evolution. The 
features we discussed are: (i) how does a possible coherence of the deconfined or hadronizing systems (modelled here 
t-H ' by some external stationary force occurring in the Langevin-like equations) describing the process of deconfinement 

■ or hadronization influence the 2-particle BEC function CiiQ = \J {p^ — P^) 2 ), where two particles are characterized 

by their four- momenta p^ and p' , and (ii) what is the true origin of the experimentally observed Q-dependence of 
the C2(Q) in the approach used in [2,1]? 

The only physical meaning of the memory term K(p^), the noise spectral function ip(p^), and the so-called coherence 
function a = a(m 2 ,p 2 ,/3) related to them [1] 
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were not yet investigated carefully (here ui is the Fourier transformed component of 4-vector p^ and it is conjugated 
to time t). We have already established [1] that at p^ — p 1 ^ ^ there is a finite volume tl(r) of a two-particle 
• • . emitter source, and the correlation picture does exist. The space in which the massive fields are settled down has 
an own characteristic length L c h ~ M^h m the same sense like the massive field with the mass m c for the Compton 
length A c = m^ 1 is served. The appearance of the 5-like distribution p(u>,e) = J dxe( %u> ~ e ') x — > 5{x) as e — > is the 
consequence of the infinite volume of produced particles (formally, e is given by the physical spectrum). The volume 
Cl(r) is conjugate to some characteristic mass scale M c h, Q(r) ~ 1/M^ h in such a way that M c /, — * oo when the 
correlation domain does not exist. Hence, M c h — > obeys the condition a — > that means no distortion (under the 
constant force P) acting on the system of produced particles. On the other hand, any distortion which disturbs this 
system by any reason (the coherence function a > 0) leads to a finite domain of produced particles (M c h 0). Rather 
a strong distortion, a — > oo, satisfies the point-like region of the particle source emitter. 

It is a difficult task to derive both K(pjj,) and V'(Pm) usm g the general properties of the quantum picture. The fact 
which is worthwhile to stress is that the knowledge of one of these functions leads to understanding of the other one, 
because of the normalization condition [2] 
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which is nothing else but the generalized fluctuation-dissipation theorem. The spectral properties 
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follow from (2) with 
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In this paper we would like to focus our attention on the role of the hadron mass influencing correlations between 
particles. To solve this problem, one has to derive the memory term K(p /1 ) using the general properties of QFT. 

It is supposed that we work with the fields corresponding to the thermal field with the standard definition of 
the Fourier transformed propagator F[G(p)] 

F[G(p)\ = G(x -y)=Tr {T[$(x)<%)] P/3 } , (5) 

with p/3 = e~P H /Tre~P H being the density matrix of a local system in the equilibrium with the temperature T — 
under the Hamiltonian H . 

We consider the interaction of <I>(x) with the external scalar field which is given by the potential U. This potential 
contrary to an electromagnetic field is a scalar one, but not a component of the four-vector. The Lagrangian density 
looks like 

L{x) = d^*(x)d^<5>{x) - (to 2 + U)<f>*{x)$(x) (6) 

and the equation of motion is 

(V 2 +m 2 )$(ar) = -J(x), (7) 

where J(x) — U$>(x) is the source density operator. Such a simple model allows one to investigate the origin of 
the occurrence of the condensate in the static restricted potential of the external field. We are interested in the 
origin of the unstable state of the thermalized equilibrium in the nonhomogencous external field under the influence 
of the source density operator J(x) — UQ(x). For example, the source can be considered as the 5- like generalized 
function J(x) — (lp(x, e)<f>(x), where p(x, e) is the 5-\ike succession giving the <5-function as e — > {jl is some massive 
parameter). This model is useful because the 5-\ike potential U(x) provides the model conditions to restrict the 
particle emission domain or the deconfinement region. We suggest the following form 

J(x) = -H(id lt )9{x) + J R {x), (8) 

where the source J(x) is decomposed into the regular systematic motion part E(ic^) $(x) and the random source 
Jr(x). Thus, the equation of motion (7) becomes 

(V 2 + to 2 - Z(id^)Mx) = -J R (x), (9) 

and the propagator satisfies the following equation: 

\- V l + m 2 - E(p M )]G(p„) = 1. (10) 

Now let us introduce the general non-Fock representation of the canonical commutation relation (CCR). To this, 
we consider the operator generalized functions 

b(x) =a(x)+R(x), (11) 

b+(x) = a+(x) +R + (x), (12) 
where the operators a(x) and a + (x) obey the CCR-relations: 

[a(x),a(x')} = [a + (x),a + (x')} = 0, (13) 

[a(x),a+(x')]=5(x-x'). (14) 

The operator generalized functions R carry random features describing an action of the external forces. 

Both b + and b obviously define the CCR representation. For each function / from the space S'(3?oo) of smooth 
decreasing functions one can establish new operators &(/) and b + (f) 

b(f) = J f(x)b(x)dx = a(f) + J f(x)R(x)dx, (15) 
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b+(f) = I f(x)b+(x)dx = a+(f) + J f(x)R+(x)dx. (16) 

The transition from the operators a(x) and a + (x) to b{x) and b + {x), obeying those commutation relations as a(x) and 
a + (x) leads to linear canonical representations. If both &(/) and b + (f) create the Fock representation of the CCR, 
then one can find the operator U obeying the following conditions: 

Uaif)!!- 1 = &(/), (17) 

Ua+ifiU- 1 = &+(/). (18) 

Now we are going to a simple physical pattern. Let us define the differential evolution (in time) equation, where 
the sharp and chaotically fluctuating function, obeying this equation, is the main object. Since we deal with the 
continuous time, this allows one to formulate the stochastic differential equation applied to each (analytical) function 
under the distortion of a random force. In classical mechanics, the stochastic processes in a dynamic system are under 
a weak action of the "large" system [7]. The "small" and the "large" systems are understood in the sense that the 
number of the states of freedom of the first system is less than that of the other one. We do not exclude even the 
interplay between these systems. In the case when the "large" system is in equilibrium state (in may be a thermostat 
state) our method will allow one to describe the approximation to the distribution of the probability to find physical 
states in the " small" system. On the quantum level the role of such a " small" system is played by the restricted region 
of produced particles, a particle source emitter with a definite size, which we study in this paper. 

Following the idea of classical Brownian motion [8,3] of a particle with a unit mass, a charge g and velocity v(t) 
in the external, let say, electric field E, one can write the formal equation describing the evolution (in time) of this 
particle in the following form: 

d t v(t) =-yv(t)+F + gE, (19) 

where F stands for a random force subject to the Gaussian white noise; 7 is a friction coefficient. 

Referring to [2,1] for details, let us recapitulate here the main points of our approach in the quantum case. The 
collision process produces a lot of particles out of which we select only one (we assume for simplicity that we deal only 
with identical bosons) and describe it by stochastic operators b{p, t) and b + (j>, t), carrying the features of annihilation 
and creation operators, respectively (the notation is the usual one: p is 3-momentum and t is a real time). The 
rest of the particles are then assumed to form a kind of heat bath, which remains in equilibrium characterized by 
a temperature T (which will be one of our parameters). We shall also allow for some external (to the above heat 
bath) influence to our system. The time evolution of such a system is then assumed to be given by a Langcvin-type 
equation [2,1] for the new stochastic operator b(p, t) 

id t b(p,t)=A(p,t) + F(p,t)+P (20) 

(and a similar conjugate equation for b + (p,t)). We assume that an asymptotic free undistorted operator is a(p, i), 
and the deviation from the asymptotic free state is provided by the random operator R(p,t): a(p,t) — > b(p, t) = 
a(p,t) + R(p,t). It means that, e.g., the number of particle density (a physical number) (n(p, t)} ph (n(p)) + O(e), 
where (n(p,t)) ph means an expectation value by a physical state, while (n(p)) denotes that by an asymptotic state. 
In case we ignore the deviation from the asymptotic state in equilibrium, one obtains the ideal fluid. Otherwise, 
one has to take into account the dissipation term, and this is the reason that we use the Langevin scheme to derive 
the evolution equation, but only on the quantum level. We derive the evolution equation in the integral form which 
reveals effects of thermalization. 

Equation (20) is supposed to model all aspects of the deconfincmcnt or hadronization processes. The combination 
A(p, t)+F(p, t) represents the so-called Langevin force and is therefore responsible for the internal dynamics of particle 
emission in the following manner: the memory term A causes dissipation and is related to stochastic dissipative forces 
[2] 

/+OC 
dTK(p,t-T)b(p,T), (21) 
-00 

with K (p, t) being the kernel operator describing the virtual transitions from one (particle) mode to another. The 
operator F{p 1 1) is responsible for the action of a heat bath of absolute temperature T on a particle in the heat bath, 
and under the appropriate circumstances is given by 
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Our heat bath is represented by an ensemble of coupled oscillators, each described by the operator c(p^) such that 
[c(p t _ l ),c + (p' fl )] 5 4 (p ll — p'^), and characterized by the noise spectral function V(jv) [2,1]- Here, the only statistical 
assumption is that the heat bath is canonically distributed. The oscillators are coupled to a particle which is in turn 
acted upon by an outside force. Finally, the constant term P in (20) (representing an external source term in the 
Langevin equation) denotes a possible influence of some external force. This force would result, for example, in a 
strong ordering of phases leading, therefore, to a coherence effect. 
The solution of equation (20) is given in by 

Kp„) = \t7— IHp») + P(«»P, <0], (23) 

where u> in p(u, e) was replaced by the new scale uop = lu/P. It should be stressed that the term containing p(tup, e) 
as e — > gives the general solution to equation (20). Notice that the distribution p(uip,e) indicates the continuous 
character of the spectrum, and the arbitrary small quantity e can be defined by the special physical conditions or the 
physical spectra. On the other hand, this p(ujp,e) can be understand as the temperature-dependent succession (4) 
where e — ► /3 _1 . Such a succession gives the restriction on the /3-depcndent second term in the solution (23) where at 
small enough T there will be a narrow peak at w = 0. 

^From the scattering matrix point of view the solution (23) has the following physical meaning: at enough outgoing 
past and future, the fields described by the operators a(p^) are free and, thus, the initial and the final states of 
dynamic system are both characterized by the constant amplitudes of states. Both of these states <p(— oo) and ip(+oo) 
are related to each other by some operator S(R) carrying out the transformation of the state <p(— oo) to the state 
<p(+oo) and depending on the behaviour of 

^(+oo) = <p(R) = S(R)<p(-oo). 

In accordance with this definition, it is natural to identify S(R) as the scattering matrix in the case of arbitrary 
sources giving rise to the intensity of R. 

Based on QFT point of view relation (11) indicates the appearance of the terms containing nonquantum fields 
which are characterized by the operators R(p^). Hence, there will be terms with R in the matrix elements, and these 
R cannot be realized via real particles. The operator function R{Pfi) could be considered as a limit of an average value 
of some quantum operator (or even a set of operators) with intcsity increasing up to infinity. The later statement 
could be realized in the following mathematical representation [2]: 



Rip,,) = yJaElp^pp), 



where a is the coherence function which gives the strength of the average 'B(p ll ,p ll ). We shall find this coherence 
function carrying the dependence on (3, the particle mass m and the chemical potential p. 

In principal, the interaction with the fields described by R is provided by the virtual particles, the process of 
propagating of which is given by the potentials defined by the R operator function. 

The condition M c h — » (or fl(r) ~ > oo) in the representation 



lim "Eip^p' ) Jim ft(r) n(u, (3) exp(-<7 2 /2) -> -L- n(u,0), 
v^K Q 2 ->0 M ch 



with [2] 



means that the role of the arbitrary source characterized by the operator function R{p^) in b(p fl ) = a{p^) + R{p ti ) 
disappears. 

Let us go to the thermal field operator $(x) by means of the linear combination of the frequency parts (f> + (x) and 

*{x) = ±y> + (x)+4>-{xj\ (24) 
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composed of the operators and 6 + (p M ) as the solutions of equation (20) and conjugate to it, respectively: 

d 3 p 



cj) (x) 



(2 7 r) 3 2(p 2 + m 2 ) 1 / 



(25) 



d 3 P 



(26) 



One can easily find two equations of motion for the Fourier transformed operators b{p ll ) and b + {p fl ) in S^R^) 

[lu - K( Pft )]b(p^) = F( Pft ) + p(u P , e), (27) 



[w - X + (^)]6+(^)F + ( Pai ) + e), (28) 

which are transformed into new equations for the frequency parts <j) + {x) and 4>~(x) of the field operator <&{x) (24) 



(29) 



-i^-(a;)- / /f+(a;- ? /) ( / ) -(y)rfy J F 1 +( a ;)+P9o J D(a ; )e" et . 
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(30) 



Here, the field components 4> + {x) and _ (a;) are nonlocalized under the effect of the invariant formfactors K{x—y) and 
K + (x — y), respectively. In general, these formfactors can admit the description of locality for nonlocal interactions. 
The function D(x) in (29) and (30) obeys the commutation relation 



[*(a;), $(»)]_ = -iD(x) 



and looks like [9] 



D(x) = — e(x°) 

Z IT 



S(x 2 ) 



m 



2 J 'x* 



(31) 



(32) 



where e(x°) and 6(x 2 ) are the standard unit and the step functions, respectively, while J\{x) is the Bcsscl function. 
On the mass-shell D(x) becomes 



D(x)*—e(x°) 

Z 7T 



«5(x 2 )-^6(x 2 ) 



(33) 



At this stage, it is needed to stress that we have got new generalized evolution equations (29) and (30) keeping the 
general features of propagating and interacting of the quantum fields with mass m in the heat bath (reservoir) and 
chaotically distorted by other fields. To proceed to a further analysis, let us rewrite the system of equations (29) and 
(30) in the following form: 



ido<f> + (x) - K(x)*<j>+(x) = f(x), 



(34) 



-%d Q 4r{x)-K + {x)*4r{x) = f + {x), 

where A(x) * B(x) means the convoluted function of the generalized functions A(x) and B(x), and 

f(x)=F(x)+Pd D(x)e- £t . 



(35) 



(36) 



Applying the direct Fourier transformation to both sides of eqs. (34) and (35) with the following properties of the 
Fourier transformation: 



F[K(x)*<f> + (x)] = F[K(x)]F[<f> + (x)], 



(37) 
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one can get two equations 



[-p°-K(p„)]4>+(p tl ) = F[f(x)}, 



(38) 



\p°-K+(p„)}4>-(p,) = F[f+(x)]. 
Multiplying cqs. (38) and (39) by p° — K + (p^) and — p° — K (p M ), respectively, one finds 

[-p° K(p,)][p° K+(p,)Mp,) = T{p,), 

where 

TGv) = [p° K+(p,)]F[f(x)] [p° + K( P ,)]F[f + (x)}. 



(39) 



(40) 



(41) 



Now we are at the stage of the main strategy: one should identify the field $(x) introduced in (5) and the field $(x) 
(24) built up of the fields <fi + and <p~ as the solutions of generalized equations (29) and (30). The next step is our 
requirement that the Green's function G(p M ) in (10) and the function T(p fl ,f3), satisfying cq.(40) 



[p a + k{p^][-p + k+{p^)]f{p^) = i, 

have to be equal to each other, where [9] 



G{p^^G(p 2 l9 \m 2 ) 



■i ^ _ i-g 2 £(p 2 ,m 2 ) 

m 2 — p 2 — it 



(42) 



(43) 



with g being the scalar coupling constant and the one-loop correction of the scalar field £ << 1/m 2 at 1/4 < (ra 2 /p 2 ) < 
1. It means that we define the operator kernel K(p^) in (27) from the condition of the nonlocal coincidence of the 
Green's function G(p M ) in (10) and the thermodynamic function r(p p ,/3) from (42) in 

lirn MG(^)-f(^,/3)]=0. 
We can easily derive the kernel operator K{p^) in the case of its real nature, i.e., 



K(p fl ) = (m 2 +f)> 



where 



and 



?(p 2 ,m 2 ) = 



l + . 9 2 £(p 2 ,m 2 ) 1- 



m 2 + p 2 



2tt 



96tt 2 to 2 Vv^ 



1 ] , p 2 ~ m 2 , 



£(p 2 ,m 2 ) = 



1 
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Ml 
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V3, 



, p 2 ~ 4m 2 . 



The ultraviolet behaviour at \p 2 \ » m 2 leads to 



£(p 2 ,™ 2 ) 



-1 



32 7T 2 p 2 



lnJ£l _^_^ 6 (p 2 ) 
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(44) 



(45) 



(46) 



(47) 



Out of many details (for which we refer to [2,1]) important in our case is the fact that the 2-particle BEC function 
for like-charge particles is defined as 



C 2 (Q) = X(N) ■ -/^'-P,. = X(N) ■ [1 + D{p^)} , 
f(Pn) ■ f(P'u) 



(48) 



G 



where f(p^,p'^) = (& + (p^)b + "(p)J&(£V)Kp^)) and f(Pn) = {b + (Pn)HPn)) are trie corresponding thermal statistical 
averages with b(pff) being the corresponding Fourier transformed solution (23). The multiplicity N depending factor 
is equal to x(-^0 = {N{N — 1))/(N) 2 . As shown in [2], (notice that operators R{p^) by definition commute with 
themselves and with any other operator considered here): 

f(p»,p'„) = /(iV) • M) + (~a + (p,ya(p',)}(a + (p'^(p,)) + (a+ (pjatfj) R + + 

+ {a + {p^a{p,))R + { P ,)R{p^, (49) 

/Gv) = (« + + l^(^)l 2 - (50) 

This defines D(p IJi ,p'^) — f(p tJ .,p' l j)/[f{pn) • /(p^)] — 1 in (48) in terms of the operators a(p M ) and R{p^) which in our 
case are equal to 

5(P,)-^ and fl( P ,)-^T. (51) 

This means, therefore, that the correlation function C2{Q), as defined by eq. (48), is essentially given in terms of 
p(u)p, e) and the following two thermal averages for the thermostat operators F(p, t): 

(F+(p,t)Ftf,t'))5 3 (p-i?) J g m P ,)\ 2 nH e +,u ' t - t '», (52) 

(F(p,t)F+tf,t'))S 3 (p-i?) J g |^)| 2 [l + nH]e— <*-*'), (53) 

where n(u) — {exp [(u> — p)(i] — 1} 1 is the number of (by assumption - only bosonic in our case) oscillators of energy 
w in the reservoir characterized by parameters fi (chemical potential) and inverse temperature /3 [10]. Notice that 
with only delta functions present in (52) and (53) one would have a situation in which our deconfined or hadronizing 
system would be described by some kind of white noise only. The integrals multiplying these delta functions and 
depending on (a) momentum characteristic of a heat bath V(Pm) an d (b) assumed bosonic statistics of produced 
secondaries resulting in factors n(ui) and 1 + n(oj), respectively, bring the description of the system considered here 
closer to reality. 

It is easy to realize now that the existence of BEC, i.e., the fact that (^(Q) > 1, is strictly connected with nonzero 
values of the thermal averages (52) and (53). However, in the form presented there, they differ from zero only at one 
point, namely for Q — (i.e., for p^ — p' ). Actually, this is the price one pays for the QFT assumptions tacitly made 
here, namely for the infinite spatial extension and for the uniformity of our reservoir. However, we know from the 
experiments, e.g., [11,4-6] that C^Q) reaches its maximum at Q = and falls down towards its asymptotic value of 
Ci = 1 at large of Q (actually already at Q <~ 1 GcV/c). To reproduce the same behaviour by means of our approach, 
one has to replace the delta functions in eqs. (52) and (53) by functions with supports larger than limited to a one 
point only. This means that these functions should not be infinite at = p^ — p'^ = but remain more or less 
sharply peaked at this point, otherwise remaining finite and falling to zero at small but finite values of \Qn\ (actually 
the same as those at which C^Q) reaches unity) 

S(P, • exp[-(^ - ppL^(r)(p, - p' v )\. (54) 

Here we replaced the 5- function with the smearing (smooth) dimcnsionless generalized function Cl(q = Qr) = 
exp[— (pn — p' l f)L^ v (r)(p l/ — p' v )], where the tensor L^ u is defined by the geometry of the space and f2o has the 
same dimension as the 6 - function (actually, it is nothing else but 4-dimensional volume restricting the space-time 
region of particle production). 

In this way we tacitly introduce a new parameter, r M , a 4- vector such that it has dimension of length. This defines 
the region of nonvanishing particle density with the space-time extension of the particle emission source. Expression 
(54) has to be understood in the sense that Cl(Qr) is a function which in the limit of r — > oo becomes strictly a, 5 - 
function. With such a replacement one now has 



(l + a)(l + a') 
where 



(55) 
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Cl(q) = 7 • fi(g), 7 



n 2 (w) 
n(w)n(u/) 



(56) 



The coherence function a is another very important one which summarizes our knowledge of other than space-time 
characteristics of the particle emission source. Notice that a > only when P^O. Actually, for a = one has 



KC 2 (Q)< X (N)[l + 1 n(Qr)], 



(57) 



i.e., it is contained between the limits corresponding to very large (lower limit) and very small (upper limit) values 
of P. Because of this a plays the role of the coherence parameter. Neglecting the energy-momentum dependence of 
a and assuming that a' — a one gets the expression 



C 2 (Q) = X (N) 



1 



2a 



1 



(1 + a) 2 



Cl(q) 



(58) 



In fact, since in general a a' (due to the fact that to ^ a/ and, therefore, the number of states identified here with 
the number of particles with given energy n(w) is also different) one should rather use the general form (48) for C 2 
with details given by (55) and (56), and with a depending on the hadron mass and on such characteristics of the 
emission process as temperature T and chemical potential p occurring in definition of n{uS). Notice that eq. (58) 
differs from the usual empirical parameterization of C 2 (Q) [12,4-6], 

C 2 (Q) ~ [1 + A ■ fi(Qr)][l +aQ + ...], (59) 

which is nothing else but the Goldhaber parameterization [13] at a = 

C 2 (Q) = l + \e-® 2r2 

with < A < 1 being a free parameter adjusting the observed value of C 2 (Q = 0) which is customary called "coherence 
strength factor" or "chaoticity" meaning A = for fully coherent and A = 1 for fully incoherent sources; a is a c- 
number, and with tt(Qr) represented usually as Gaussian. Recently eq. (59) has found strong theoretical support 
expressed in detail in [14]. Coming back to the hadron mass dependence of C 2 (Q), in particular, the correlation 
radius, we find that this dependence comes from a-coherence function (1) containing the operator kernel K{p^) 
defined correctly up to the second-order of the scalar coupling constant g 2 in (44) in the framework of QFT. 



The a - representation in (1) needs to be clarified. In fact, M* h may be decomposed into two parts: M* h — > ■ 

^ is the small massive scale characterized by the time-like scale r°, while M c h being the characteristic 

l 



M c \, where M ch 



mass associated with the space inverse components r z , r t , M^ h <~ (r z r 2 ) 1 . Taking into account the properties of the 
distribution [p(wp, e)] 2 one can suggest the following replacement M^ h [p{ojp, e)] 2 — > [M^ p(0, e)][M?7i p(wp, e)], where 
the first multiplier is of the order of 0(1) while the second one reflects the massive scale of the particles with the mass 
to, i.e., M^ h p(ujp,e) - 0(m 2 ). Thus, 



a- O 



(60) 



Let us return to the problem of Q-dependence of BEC. One more remark is in order here. The problem with the 
<5(p M — p'^) function encountered in two particle distributions does not exist in the single particle distributions which 
arc in our case given by eq. (50) and which can be written as 

/(Pm) = (! + «)■ S(p M ,p M ), (61) 
where E(p^,p^) is the one-particle distribution function for the "free" (undistorted) operator a(p^), namely 



Z(pn,Pv) = (a + (p AJ )a(p M )) = n a 



uj - K{p ll ) 



n(uj). 



(62) 



Notice that the actual shape of f{p li ) is dictated by both n(ui) (calculated for fixed temperature T and chemical 
potential p at energy to as given by the Fourier transform of field operator K (44) and the shape of the reservoir in the 
momentum space provided by ip(p^)) and by the 5 - like distribution of external force p(u)p, e). On the other hand, 
it is clear from (61) that (N) = (N c h) + (N co h) (where (N c h) and {N co h) denote multiplicities of particles produced 
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chaotically and coherently, respectively. For illustration we plotted in Fig. 1 in detail (with the correct Gaussian 
shape of £l(q) function) the dependence of C2(Q) on different values of a = 0, 0.25, 1, 4 and compared it with the 
case when the second term in eq. (58) is neglected, as is the case of the majority of phenomenological fits to data. 
Finally, we present the lower bound of the particle source size 



r 2 > 



1 



■In 



7X(A0 



Q2 (l + a)(l + a')[C 2 (0)- X (N)} 



(63) 



carrying the dependence of: 

- the absolute temperature of a heat bath and the chemical potential (7- factor defined in (56)), 

- the mean multiplicity factor x(iV) defined within formula (48)), 

- the maximal value of the BEC function C 2 (Q — 0) at the origin under the following condition: 



C 2 (0) < X (N) 



(l + a)(l + a') 



(64) 



- and the hadron mass m given by a and a' defined in approximation (60) with the kernel operator K(P fl ) (44) 
calculated in QFT. 

To summarize: using the QFT supplemented by Langevin-like evolution equations (29) and (30) to describe de- 
confincment or hadronization processes we have derived the two-particle BEC function in the form explicitly showing 
the origin of both the so-called coherence (and how it influences the structure of BEC) and the Q-dependence of 
BEC represented by the correlation function C 2 (Q). The dynamic source of coherence is identified in our case with 
the existence of a constant external term P in the evolution equation. Therefore, for P — > 00 we have all phases 
aligned in the same way and C 2 (Q) = 1. This is because the coherence has already been introduced on the level of a 
particle production source as a property of fields or operators describing produced particles. It is therefore up to the 
experiment to decide which proposition is followed by nature: the simpler formula (59) or rather the more involved 
(48) together with (55). ^From Fig. 1 one can see that differences between both forms are clearly visible, especially 
for larger values of the coherence function a. 




Q [GeV] 

FIG. 1. Shapes of C2(Q) as given by eq. (58) - upper panel and for the truncated version of (58) (without the middle term) 
lower panel. Gaussian shape of 57(g) was used in both the cases. 



^From our approach it is also clear that the form of C 2 reflects distributions of the space-time separation between 
the two observed particles. 
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Finally, we would like to stress that our discussion is so far limited to only a single type of secondaries being 
produced. It is also aimed at a description of deconfincmcnt or hadronization understood as kinetic freeze-out in some 
more detailed approaches. This is enough to obtain our general goals, i.e., to explain the possible dynamic origin of 
coherence in BEC, the origin of the specific shape of the correlation C^iQ) functions and to explain the dependence 
of the correlation radius on the hadron mass which is carryed out by the coherence function a, as seen from the QFT 
perspective. It is then plausible that in the general description of the BEC effect they should be somehow combined, 
especially if experimental data indicate such a necessity. 

Part of this work is based on a collaboration with G. Wilk and O. Utuyzh. I have greatly benefited from our many 
animated discussions. I am also grateful to S. Tokar for his kind hospitality and for many fruitful discussions. 
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